Hybridization and interference effects for localized superconducting states in strong 

magnetic field 
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Within the Ginzburg - Landau model we study the critical field and temperature enhancement 
for crossing superconducting channels formed either along the sample edges or domain walls in thin- 
film magnetically coupled superconducting - ferromagnetic bilayers. The corresponding Cooper pair 
wave function can be viewed as a hybridization of two order parameter (OP) modes propagating 
along the boundaries and/or domain walls. Different momenta of hybridized OP modes result in the 
formation of vortex chains outgoing from the crossing point of these channels. Near this crossing 
point the wave functions of the modes merge giving rise to the increase in the critical temperature 
for a localized superconducting state. The origin of this critical temperature enhancement caused 
by the wave function squeezing is illustrated for a limiting case of approaching parallel boundaries 
and/or domain walls. Using both the variational method and numerical simulations we have studied 
the critical temperature dependence and OP structure vs the applied magnetic field and the angle 
between the crossing channels. 

PACS numbers: 74.25. Dw 74.78.Na 74.78.Ch 



I. INTRODUCTION 

Recent experimental and theoretical studies of ferro- 
magnet/superconductor (F/S) heterostructures have re- 
vealed a rich physics of magnetic and transport properties 
of these systems (see, e.g., Refs. [HHI for review). A con- 
siderable amount of attention in these studies has been 
devoted to the effect of ferromagnetic domain structure 
on the critical temperature of superconductivity nucle- 
ation (see, e.g., Refs. and references therein). This 
effect originates from both the exchange and electromag- 
netic (orbital) mechanisms of interaction between super- 
conducting and magnetic ordcrings. The electromagnetic 
mechanism appears to be a dominant one for the experi- 
mental works which investigate F/S bilayers designed to 
suppress the proximity effect (see, e.g., Refs. [t|). The 
non-uniform magnetic stray field of domain walls can re- 
sult in the formation of a localized Cooper pair wave 
function at temperatures exceeding the superconducting 
critical temperature in the bulk. Such localized order 
parameter (OP) distributions and corresponding phase 
diagrams have been studied in Refs. for a generic 

case of a step- like profile of a stray magnetic field in a su- 
perconducting thin film. At a fixed temperature T for a 
certain critical amplitude Bo \H C 2(T) < Bq < H c3 (T)j of 
a step - like magnetic stray field profile the structure of 
localized OP wave function coincides with the one which 
formes in the vicinity of superconductor/vacuum or su- 
perconductor/insulator flat interface in a parallel field 
equal to B , where H c2 (T) and H c3 (T) = 1.695if c2 (T) 
are the critical fields of bulk superconductivity and sur- 
face superconductivity, respectively (Toj . 

Certainly, in a restricted sample geometry the localized 
OP wave function should appear not only near the do- 
main walls but also at the film edges. Thus, for temper- 
ature values T between T c2 {B ) and T c3 (B ) > T c2 (B ), 
defined by the conditions H C 2(T C 2) = Bq and H C 3(T c3 ) = 



Bq, respectively, there appears a set of localized super- 
conducting modes propagating along the sample edges 
and/or domain walls and decaying in the transverse di- 
rection. Hereafter these localized modes will be called 
superconducting channels. One can pose a natural ques- 
tion whether the superconducting critical temperature 
T c could be increased due to the overlapping of the OP 
modes localized near the various types of superconduct- 
ing channels. Indeed, such increase in the T c value is well 
known to occur for two parallel superconductor/vacuum 
interfaces (i.e. for a superconducting film) placed in a 
homogeneous magnetic field H parallel to the interfaces 
(see, e.g., Refs. [TflHH). The critical temperature in this 
system monotonously increases with the film thickness 
decreasing and saturates at the T c0 = T C (H = 0) value in 
the limit of small thicknesses. Taking non-parallel bound- 
aries we get a superconducting wedge placed in a mag- 
netic field. The squeezing of the superconducting wave 
function at a small wedge angle is known to cause a strong 
increase in the critical temperature T c (see Refs. [l2l - [2ll ] ) . 

It is the goal of our work to analyze the effect of such 
wave function squeezing on the critical temperature en- 
hancement for various types of crossing superconducting 
channels localized near the sample edges and/or domain 
walls. Within the linearized Ginzburg - Landau model 
we carried out the calculations of the critical temperature 
of localized superconductivity for a simple hybrid system 
consisting of a thin superconducting film placed in a non- 
uniform stray magnetic field of straight domain bound- 
aries in a ferromagnetic substrate underneath. An origin 
of the critical temperature enhancement can be clearly 
demonstrated for a generic example of parallel bound- 
aries. Generalizing the textbook solution for the critical 
field of a thin superconducting film [kJ [LH we consider 
two model problems: (i) parallel domain boundary and 
sample boundary separated by the distance D; (ii) two 
parallel domain boundaries at the separation D. In both 
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cases the dependence of critical temperature T c vs the 
distance D reveals a maximum at the value of the order 
of the magnetic length Lb = y/hc/\e\Bo, where h is a 
Planck constant, e = — |e| is an electron charge. At this 
distance the wave functions of the OP modes localized 
near neighboring domain walls or near the domain wall 
and the sample edge merge and form a single supercon- 
ducting channel. Turning now to the case of non-parallel 
domain walls (which cross at the angle ipo <C t/2) or 
a domain wall crossing a sample boundary at a certain 
small angle ipg one can expect that the most favorable 
conditions for superconductivity nucleation realize at the 
distance ~ Lb^Po from the crossing point. As a result, 
for rather small angles ipo <C tt/2 the center of the en- 
ergetically favorable OP distribution is shifted from the 
crossing point. This shift is accompanied by the striking 
phenomenon: the critical temperature T c of the local- 
ized superconductivity monotonously increases when the 
tilting angle (fo tends to zero. Such behavior makes the 
dependence T c (tp ) nonanalytic at ipo = where the crit- 
ical temperature exhibits a jump to the value T C 3 (T^) 
for the case of the domain wall crossing the sample edge 
(for two crossing domain walls) . 

In our work we suggest a simple variational procedure 
which allows us to get an approximate solution of a lin- 
earized Ginzburg - Landau equation describing the hy- 
bridization of the localized superconducting states for an 
arbitrary tpo angle. We hnd out that different momenta 
of hybridized OP modes arc responsible for the formation 
of vortex chains outgoing from the crossing point of the 
channels. The effect of these chains on the OP trial func- 
tion appears to be extremely important for rather large 
ipo angles close to tt/2. The important effect of such 
vortex chain can be illustrated for the simplest example 
of the superconductivity nucleation at the wedge corner. 
Our calculations show that only taking into account the 
vortex chain one can get a proper crossover to the H c s 
field with the increase in the wedge angle up to the fiat 
one. The solution of a generic problem describing the su- 
perconductivity nucleation in a wedge allows us to find 
appropriate trial wave functions for superconducting OP 
nucleating near the domain walls intersecting the film 
edge or near the crossing point of domain boundaries. 
The change in the <po angle is found to be accompanied 
by the change in the orientation of vortex chains and the 
intervortex distance. Our analytical findings based on 
the variational procedure are in a good agreement with 
direct numerical simulations. 

Applying an external magnetic field perpendicular to 
the film plane one can observe the increase in the critical 
temperature of the domain wall superconductivity due 
to the partial magnetic field compensation inside the do- 
mains. Using both variational and numerical approaches 
we have composed the phase diagram of the F/S bilayer 
in the plane temperature T - external magnetic field H 
and discuss the transitions between different supercon- 
ducting states by varying temperature T, external mag- 
netic field H and magnetic stray field amplitude Bq. 



The paper is organized as follows. In Sec. II we mainly 
focus on the variational analysis of the superconductivity 
nucleation and interference of localized superconducting 
states. In Sec. Ill we present the results of direct numer- 
ical simulations which support our analytical findings. 
Finally, the results are summarized in Sec. IV. 



II. VARIATIONAL APPROACH 

We start our study of the hybridization and inter- 
ference effects for interacting superconducting channels 
forming in a thin- film magnetically coupled F/S bilayer 
with the GL variational procedure focusing on the anal- 
ysis of appropriate trial wave functions. Let us analyze 
the problem of the OP nucleation in a thin supercon- 
ducting film placed in the non-uniform magnetic field 
B(r) = Hz n + b(r) induced by the external sources and 
the magnetic domain walls in a ferromagnetic substrate, 
respectively. Note that we will restrict ourselves to the 
case of a step-like distribution of the magnetic stray fields 
of the domain walls and neglect the effect of the magnetic 
field components parallel to the film plane. In particu- 
lar, for a single domain wall we take b(r) = zoi?o s ig n (a?)i 
where the x axis is directed perpendicular to the corre- 
sponding domain wall. Thus, we neglect the deviations 
from this step-like field model caused by a finite thick- 
ness of a superconducting film and decay of the magnetic 
stray field at large distances from the domain walls (see 
discussion in Refs. 1, [H H3). 



A. Linearized Ginzburg — Landau model 

The superconducting critical temperature T c can be 
routinely determined from the linearized GL equation 
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as the highest possible value: T c = max T, corresponding 
to the lowest "energy level" E oc £ _2 (T) of the eigenvalue 
problem (fTJ). Here ^(r) is the OP distribution, A(r) is 
the vector potential corresponding to the total magnetic 
field B(r) , m is the electron mass, £ (T) = £o/ ^/l — T/T c q 
is the superconducting coherence length. Alternatively, 
T c is known to be determined from the variational prob- 
lem 
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and the integration is performed over the superconduct- 
ing volume. The wave function ^(v) satisfies the bound- 
ary condition 
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where n is a unit vector normal to the sample boundary 

r. 

The complex- valued wave function ^E'(r) can be written 
in the form \I/(r) = /(r)e ie( - r \ where /(r) and O(r) are 
the absolute value and the phase of 'F(r), respectively. 
Thus, we can rewrite Eq. ([2]) in the form 
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j {4m 2 v 2 (r)/ 2 (r) + ft 2 [V/(r)] 2 } 
h 2 J f 2 (r)d 2 r 



d 2 r 



(4) 



where v s (r) = [W8(r) — 2eA(r)/c]/2m is a superfluid 
velocity. 



B. Localized superconducting modes forming in 
the presence of parallel boundaries 

We begin with the consideration of the superconduc- 
tivity nucleation for parallel sample edges and/or do- 
main walls in a thin- film F/S bilayer. Let us choose 
the x (y) axis to be perpendicular (parallel) to these 
boundaries and take the gauge A = A y (x) yo . The 
Schrodinger-like equation ([!} docs not depend on y coor- 
dinate and one can generally find the solution in the form 
ty{x,y) = fic(x)exp(—iky), where the function fk{x) 
should be determined from the following eigenvalue prob- 
lem: 



d 2 fk(x) 
dx 2 



+ U(x)f k (x) 
U(x) = 
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We calculate the dependence of the shift of the crit- 
ical temperature (1 — T c /T c o) on the distance D using 
both numerical solution of the eigenvalue problem J5]) 
and variational approach. 

Two superconductor/vacuum interfaces. For two par- 
allel superconductor/vacuum interfaces positioned at x = 
±D/2 and forming a superconducting slab of the finite 
thickness D placed in a homogeneous parallel magnetic 
field B z = H we should apply the boundary condi- 
tion dfk(x)/dx\ x= ±D/2 = 0. The dependence of E = 
(1 - T c /T c0 ) (L H /£o) 2 on D (where L H = y/hc/\eH\ is 
a characteristic length scale in a magnetic field), origi- 
nally calculated by Saint- James and de Gennes [10|, UM 
is shown in Fig. [T] by the solid line (1). The critical tem- 
perature T c tends to T c0 at D <C Ljj and, thus, the mag- 
netic field has a negligible effect on the OP nucleation for 
rather small D values. In this limit a simple approxima- 
tion for the dependence E{D) can be found provided we 
choose Ay(x) = Hx and assume the OP wave function 
to be uniform across the slab. The minimum in the de- 
pendence E(k) in this case corresponds to k = 0. Such 
approximation while being rigorously justified only for 
D L H appears to describe the monotonous decrease 
in T c with the increasing D distance even for D ~ Lh (see 
the dashed line (1) in Fig.[T]). In the limit D 3> Lh, there 




FIG. 1: (Color online) The shift of the energy E = (1 - 
T c /T c o) x (i/Co) 2 vs the distance D, obtained numerically 
(solid lines) and using the trial function approach (dashed 
lines). The curves (1), (2), (3) correspond to three model 
problems: (i) two superconductor/vacuum interfaces , (ii) two 
domain walls, (iii) domain wall and superconductor/vacuum 
interface, respectively. Here L = Lh for the first problem, 
L — Lb for the last two problems. 



are two symmetrical minima in the dependence E(k) po- 
sitioned at nonzero k values. These two solutions cor- 
respond to the same critical temperature of surface 
superconductivity: T c ~ T c3 = T c0 (1 - 0,m\H\/H^), 

where H^f = *o/2<o is the 

upper critical field at 
T = 0, and $o = ^hc/\e\ is the magnetic flux quantum. 
The corresponding OP wave functions describe the local- 
ized superconducting modes running along the bound- 
aries x = —D/2 and x = D/2. 

Two domain walls. We proceed our consideration with 
the case of two parallel domain walls separated by the 
distance D: B z = -B for -D/2 < x < D/2 and 
B z = Bq elsewhere (here we assume H = 0). In the 
limit D = the critical temperature T c equals to the 
bulk critical temperature T C 2 = T c q (1 — Bq/H^ ) in the 
uniform magnetic field Bq. In the opposite limiting case 
D > L B = yJhc/\e\B the OP wave function localized 
at the discontinuities of the magnetic field component 
B z saturates the critical temperature T c3 . What is less 
intuitively clear is that the transition from T C 2 to T C 3 
with the increasing distance D occurs via the T c en- 
hancement. Such non- monotonous T c behavior can be 
captured with a good accuracy by minimizing the en- 
ergy functional ^ for the Gaussian-like trial function: 
fk(x) = cxp(— x 2 /a 2 ), where a is a variational parame- 
ter. The choice A y (—x) = — A y (x) automatically yields 
k = 0. Then it is clear that the rise in T c is caused by 
the increase in the width of the effective potential well 
U (x) = A 2 (x) for increasing D and by the lowering of the 
ground energy level in a wider potential well. The com- 
parison of the results of the numerical solution of Eq. §5§ 
[solid line (2)] with the trial function approach [dashed 
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line (2)] is presented in Fig. Q] 

Domain wall and superconductor/vacuum interface. 
Finally, we analyze the case of a single domain wall par- 
allel to the superconducting film edge and positioned at a 
distance D from this edge. We also assume H = so that 
the only characteristic length scale is Lb ■ Similarly to the 
previous case the T c value changes non-monotonuously 
as a function of D and in both limits (D -C Lb and 
D ^> Lb) the critical temperature tends to the critical 
temperature of surface superconductivity: T c — > T c3 . At 
D — > oo there are two independent superconducting OP 
nuclei located near the surface (x = —D/2) and at the 
domain wall (x = D/2) both characterized by the same 
critical temperature T c ^. Analogously to the previous 
case of two parallel domain walls the non-monotonous 
behavior of T c is caused by the increasing width of the 
potential well U(x) = [2nA y (x)/^ - k} 2 in the Eq. (©. 
The problem can be apparently mapped on the one in 
the infinite superconducting slab in the magnetic field 
of three parallel domain walls placed at x = ±D/2 and 
x = —ZD/2. The resulting magnetic field distribution 
is an odd function for reflection respective to the plane 
x = —D/2. For such magnetic field configuration the 
ground state solution should possess the reflection sym- 
metry: f k (x - D/2) = f k (-x - D/2). Therefore for 
D < Lb we can choose the trial function in the form 
fk{x) = exp[— (x — D /2) 2 /a 2 ], where a is a variational 
parameter. Unlike the previous case the vector potential 
now is an even function A y (x — D/2) = A y (—x — D/2) 
which means that the ground state solution corresponds 
to k 7^ 0. Minimizing the energy functional over the 
parameters a and k, we obtain the shift of the critical 
temperature shown by dashed line (3) in Fig. [TJ This 
plot again demonstrates an excellent agreement with the 
numerical result [solid line (3)] for rather small D values. 

Thus the solution of these model problems allows us 
to make an important observation about the possibility 
to get the critical temperature enhancement for a pair of 
approaching boundaries or domain walls. 



C. Hybridization of localized superconducting 
modes propagating along the channels crossing at 
small angles 

Considering the problem of superconductivity nucle- 
ation for crossing boundaries or domain walls it is natural 
to start from the simplest case of small crossing angles 
when the distance D between the crossing boundaries 
changes adiabatically. 

Two superconductor/vacuum interfaces. The solution 
for two crossing superconductor/ vacuum interfaces which 
form a superconducting wedge with small corner an- 
gle x *^ 7r /2 can be found in Refs. [l2l - |2"iT ] Let us in- 
troduce a cylindrical coordinate system (r, <p, z). The 
monotonous increase in T c with the decreasing distance 
D discussed above allows to assume that the maximum 
of the OP wave function should be positioned at the 



wedge vertex (r = 0) . Substituting the simplest isotropic 
wave function ^(r) = exp(— r 2 \eH\x/2\/3hc) in the func- 
tional ([2]) and carrying out the variational procedure one 
can find the following asymptotical expression for the 
critical mag netic field H™ 3 ~ V3H c2 /x 

suppressing the 

localized superconductivity. 

ft 




FIG. 2: (Color online) Schematic view of a F/S bilayer: a 
semi-infinite thin superconducting film and a straight domain 
wall (dashed line) oriented at a certain angle tpo <C 7r/2 with 
respect to the film edge V. Bo is a stray field amplitude inside 
the domain. 

Domain wall crossing the superconductor/vacuum in- 
terface. We proceed now with the case of a single domain 
wall oriented at a rather small angle ipo with respect to 
the superconducting film edge T (see Fig. [2]). We re- 
strict ourselves by a particular case of zero external field 
H = 0. Thus, we consider the variational problem ©, 
([3]) for a semi-infinite superconducting thin film in a mag- 
netic stray field b z (r) = i?osign(a;). We choose the trial 
function /(r) in the form: 

f(r) = e Sv 2 e -0(*-*o)\ ( 6 ) 

where 5, (3 and xq are the variational parameters. Here 
we allow the center of a superconducting nucleus to be 
shift ed from the crossing point of the domain wall and the 
sample edge T only along the x axis (xo 7^ 0). The shift 
of the nucleus center is a direct consequence of the exis- 
tence of the maximum in the dependence of the critical 
temperature T c vs the separation D between the bound- 
aries. According to the above treatment this maximum 
corresponds to the D value of the order of the magnetic 
length Lb- Thus, for a small angle ipo the peak in the 
OP wave function should appear at the distance ~ Lb^Po 
from the crossing point. 

Accounting the magnetic stray field b z (r) of the do- 
main wall we take the superfluid velocity v s (r) in each 
magnetic domain to be equal to the vector potential A(r) 
chosen in the radial gauge: 

v s (r) = B r(tp - ■&) x \e\/(mc) r , 
< f < (fo , 

v a (r) = B r(-tp + 2ip - t?) x \e\/(mc) r , 

<A) < f < 7T , (7) 
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where & is the variational parameter (0 < & < tpo). 
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H C3= 1 - 695 H C2 
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FIG. 3: (Color online) The critical amplitude Bg of the do- 
main stray held of localized superconductivity nucleation in a 
thin- him semi-inhnite F/S bilayer as a function of the titling 
angle ipo- Dash-dot lines (1) and (3) are obtained from the 
variational analysis carried out both for small and large cor- 
ner angles, respectively. The dependence Bo(ipo) found from 
numerical simulations is shown by the solid line (2). 



(a) (b) 




% (degree) % (degree) 

FIG. 4: (a) Variational parameter xq vs ipo obtained using the 
trial function approach (squares). The solid line corresponds 
to the analytical dependence xo = OAir/ipo. (b) The depen- 
dences of P X 180° (squares) and i9 (circles) vs ipo obtained 
using the trial function approach. The solid lines (1) and (2) 
correspond to the analytical dependences f3 x 180° = ipo and 
$ = ipo /2.2, respectively. 

Substituting the expressions © and into the Eq. dJ) 
and carrying out the minimization over the parameters 
<5, j3, xq, and &, we derive the dependence of the criti- 
cal field amplitude Bq vs the ipo angle for the limit of 
small angles (fo -C 1. Typical plot of this dependence is 
shown in Fig. [3] by a dash-dot line (1). For ipo <C 1 this 
plot demonstrates a good agreement with our numeri- 
cal simulations [see the solid line (2) in Fig. [5] carried 
out below. It is interesting to note that the position xq 
of the OP maximum diverges inversely proportional to 
tpo as tpo — >• in accordance with our expectation [see 
Fig. IDJa)]. The parameter <5 is found to be independent 
of ipo: 6 = 0.16£o- The dependences of the variational 
parameters (3 and $ vs the tpo angle are presented in 
Fig. lib). 

One can see that even for ipo <C 1 the critical field Bq 
do differs from the H c ^ value corresponding to conven- 
tional surface superconductivity in a homogeneous field: 




FIG. 5: (Color online) Schematic view of a F/S bilayer: a thin 
superconducting him and two domain walls (dashed lines) 
crossing at a certain angle ipo <C tt/2. Bo is a stray held 
amplitude inside the domain. 



Bq — > 2.8H C 2 at tpo — > 0. At tp = the dependence 
Bq((Po) exhibits a step-like jump to the H c ^ field. This 
jump should be, of course, smeared for a finite size sam- 
ple when xq becomes comparable to the sample size. 

It is important to note that the threshold value Bq 
allows to determine the critical temperature at H = 0: 
T c = T cQ {1-B q /B*q). 

Two domain walls. The consideration in previous sub- 
sections can be easily generalized for the case of two do- 
main walls which cross at rather small angle ipo <C 7r/2 
(see Fig. [5|). Due to the symmetry of the magnetic 
field distribution we can consider only the half-space: 
< if < 7r. We take the superfluid velocity v s (r) in 
the form ([7|) and choose the trial function /(r) like in the 
following 

/(r) = e - 5 (y-yo) 2 e -P^-^) 2 ( 8 ) 

Here 5, (3, xq and yo are the variational parameters. It is 
important to note that the center of a superconducting 
nucleus is shifted from the crossing point of two domain 
walls along both the x and y axes. Substituting the ex- 
pressions ([7]) and (8]) into the Eq. ((4]) and carrying out 
the minimization over the parameters r5, f3, xq, yo, and i?, 
we derive the dependence B^ipo) for the limit of small 
angles tpo *C 1- This critical field dependence appears to 
be very close to the one shown in Fig. [3] by the dash-dot 
line (1): the relative deviation is of the order of several 
percents. We emphasize that the nucleus center (a;o,yo) 
shifts along the bisectrix of the ipo angle. 
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D. Interference of localized superconducting modes 
propagating along the channels crossing at large 
angles 

We proceed now with the consideration of the local- 
ized superconductivity nucleation for the superconduct- 
ing channels crossing at rather large angles. 

Two superconductor/vacuum interfaces. We start with 
an exemplary problem of superconductivity nucleation in 
a wedge with the corner angle \ = 71 placed in a uniform 
magnetic field H. The superconducting wave function 
localized near the wedge vertex can be considered as an 
overlapping of the superconducting modes propagating 
along the wedge sides. These modes are characterized by 
the wave vectors parallel to the different sides and, thus, 
the interference of these localized waves should result in a 
formation of vortices at the bisectrix of the wedge angle. 
The superconductivity in a wedge appears for H > H c3 
and, thus, the OP decays with the increasing distance 
from the wedge vertex. We would like to note that the 




FIG. 6: (Color online) Cross section of a superconducting 
wedge with an arbitrary corner angle \ in a uniform mag- 
netic field H . The solid circles show the vortex chain out- 
going from the wedge vertex along its bisectrix. The lines 
of a superfluid velocity v s and the contour lines of /(r) are 
presented. Dashed lines (x = t and y = t) correspond to 
v s = 0. 

interference effect and the resulting formation of a vortex 
chain have been disregarded in previous works fl2T - [l6l . l22j 
considering different types of trial functions for the OP in 
a wedge. As a consequence, all these variational calcula- 
tions provided poor agreement with the numerical results 
(see, e.g., Ref. for rather large wedge angles up to 
7r. We will demonstrate that only accounting the vortex 
chain outgoing from the wedge vertex one can obtain a 
proper crossover to the H c3 field at \ — > 7r. 

Considering a superconducting wedge (see Fig. [6]) 
and introducing the dimensionless variables in the func- 
tional (U) we come to the following expression 

1 J{vf(r');V) + [V,/-(r')] 2 }rfV 

h Jf 2 (r')d 2 r' ' 1 ' 

where h' 1 = L 2 H /Q X fl - T/T c0 ), d 2 v' = L H 2 d 2 r, 
v(.(r') = VO(r') + a(r') is a dimensionless superfluid ve- 



locity and a(r') = |77| _1 A(r') is a dimensionless vector 
potential. 

Similarly to Ref. [lH we choose the absolute value /(r') 
of Cooper pair wave function \l/(r') in the form 

/(r') = cxp{- a / 2 [l + 7sin 2 (7T0/ X )]} , (10) 

where a and 7 are the variational parameters, r' = L~^r, 
(j) = ip — x/2. To describe the vortex chain positioned 
along the bisectrix we introduce a cut along the line = 
where the superfluid velocity experiences a discontinuity. 
The tangential jump in the v s value corresponds to a con- 
tinuously distributed vorticity along this cut. For this 
purpose we divide the wedge into the angular domains 
— x/2 _ 4> < and < (f> < x/2 and take the dimen- 
sionless superfluid velocity v' s (r') in each domain equal 
to the vector potential a(r') chosen parallel to the wedge 
sides at <fi = Tx/2, correspondingly. Both gauges of the 
vector potential a(r') correspond to a uniform field H: 

v' s (r') = (t-y')x , - X /2<0<O, 
v' s (r') = (x-t)y , O<0< x /2, (11) 

where t is a variational parameter, (.t, y) is a new refer- 
ence system rotated at the angle x ~ 71 / 2 with respect 
to the original system (x, y) in the counter-clockwise di- 
rection (see Fig. [5]). Such choice of the superfluid veloc- 
ity provides a correct asymptotical behavior of the OP 
modes propagating along different wedge sides at large 
distances from the vertex. Of course, a single-valued 
wave function should vanish at the cut and, thus, the 
choice of superfluid velocity in the form Eq. (|llj) is not 
adequate for description of the wave function behavior 
close to the wedge vertex (r < £) where the absolute 
value of the OP is essentially nonzero for all angles </>. 
As a consequence, we can use this method only for wave 
functions strongly elongated along the wedge sides when 
the region close to the vertex provides a small contribu- 
tion to the functional ((9|). We will see below that this 
condition appears to break down for small angles x when 
the wave function is almost isotropic for all distances r. 

Rewriting the expressions (jTTJ) in terms of r' and 4> and 
substituting them into the Eq. © one obtains the func- 
tion hr 1 = h~ l (a, 7, t, x)- Carrying out the minimiza- 
tion over the variational parameters a, 7 and t we find 
the critical field HJ^T) = H c2 {T) x [min a)7)t (l//i)]- 1 of 
superconductivity nucleation for different wedge angles 
X- Typical plot of the dependence H™ 3 (x) is presented in 
Fig. [7] by the solid line. 

One can see that for x = n we obtain the result 
found previously by Saint- James and dc Gennes [Ioj |: 
H™ 3 (ir) ~ 1.695-ff C 2. For a particular case x ~ 7r /2 
we find H™ 3 (ir/2) ~ 2H C 2- This value appears to be in 
a good agreement with the numerical calculations car- 
ried out in Ref. [l8| (see also the comparison with nu- 
merical simulations below). The variational parameters 
corresponding to this case should be taken as follows: 
a(7r/2) = 0.14, j(n/2) = -0.5, i(7r/2) = 0.6. 

Analyzing the typical contour plots of the f(r') func- 
tion in the insets of Fig. [7] for different x on e can see that 
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FIG. 7: (Color online) (a) Typical plots of the nucleation field 
H™ 3 in a wedge vs the corner angle x- The solid line shows 
the results of our variational calculations while the dashed line 
corresponds to asymptotic behavior of the nucleation field H™ 3 
at small corner angles according to Refs. (b) 

The contour plots of f(r') function inside the areas r' < 5 for 
different x angles: O.Olvr (A), 0.5tt (B) and 0.9tt (C). 



the angular anisotropy of the OP vanishes (7 0) in 
the limit x ~ * 0. For such small angles x we find out: 
H™ 3 /H C 2 ~ 1.56v3x _1 - This asymptotical behavior de- 
viates from the correct dependence H™ 3 /H c2 = v^X -1 
found previously in Refs. [I2I - I2TI ]. This deviation for 
small angles is a natural consequence of the wave function 
isotropy as it is discussed above. With the increase in the 
wedge angle x the anisotropy parameter I7I grows and the 
wave function becomes elongated along the wedge sides 
which restores the validity of our approach. The vortex 
- free trial functions considered in Refs. (T2I - TT6L [22j can 
no more provide a correct behavior of the upper critical 
field H^. 

Domain wall crossing the superconductor /vacuum in- 
terface. The variational procedure described in the lat- 
ter subsection can be easily generalized for other model 
systems containing crossing superconducting channels. 
We proceed now with the consideration of the localized 
superconductivity nucleation in a thin-film semi-infinite 
F/S bilayer with a domain wall crossing the sample edge 
T at rather large angle tpo < it/ 2 (see Fig. [5|. To gen- 
eralize the considered approach one needs to introduce 
two vortex chains outgoing from the crossing point of a 
domain wall and the sample edge T. It is naturally to 
expect these vortex chains to be oriented along the bi- 
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FIG. 8: (Color online) Schematic view of a F/S bilayer: a 
semi-infinite thin superconducting film and a straight domain 
wall (dashed line) oriented at a certain angle <^o with respect 
to the film edge Y. Bq is a stray field amplitude inside the 
domain, H is an external magnetic field, the x axis is chosen 
to be perpendicular to the domain wall. The vortex chain is 
shown by the solid circles while the opened circles indicate 
the antivortices. 



sectrices of two wedge-shaped regions located to the left 
ant to the right of the domain wall. We follow the proce- 
dure described in the previous subsection and replace the 
vortex chains by the cuts for superfluid velocity, i.e., by 
the lines with a continuously distributed vorticity. The 
absolute value of the OP wave function and the dimen- 
sionless superfluid velocity are taken in the forms: 

/i(r) = cxp{-air 2 [l + 71 cos 2 (tt</j/¥>o)]} , 

0<Lp<(p Q , 
/ 2 (r) = 



cxp I — a 2 r 

v s (r) = 
v s (r) = 
v s (r) = 



+ 72 cos 2 (tt(<p - ip )/(ir - tpof) j, 

<A) < V < * ,(12) 



(H + B a )(r) — y) x \e\/{mc)x , 
< V < W 2 : 
(H + B )(x-r)) x |e|/(mc)y ■ 

(po/2 <(fi<Lp . 

(H - B )(x + r)) x |e|/(mc)y ■ 
</?o < if < tt/2 + (fia/2 . 
(H - B ){ri -y)x |e|/(mc)x . 

7r/2 + ipo/2 < tp <ir 



(13) 



where a±, 71, 72, and 77 are the variational parame- 
ters. The continuity condition for the wave function 'F(r) 
at ip = tpo gives the relation 



+ 7i) = 012(1 + 72) 



(14) 



For simplicity wc start with a particular case of the do- 
main wall perpendicular to the film edge T (ipo = t/2). 
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FIG. 9: (Color online) (a) Typical plots of the critical temperature T c * 3 (i/) for localized superconductivity in a thin-film semi- 
infinite F/S bilayer with a domain wall oriented perpendicular to the sample edge F (ipo = ft/2, see Fig. The solid line 
shows the results of our variational calculations while the dashed line corresponds to our numerical simulations carried out 
for a rectangular superconducting film with lateral dimensions 20£o x 20£o- The dotted lines correspond to the dependence 
of the shifted upper critical field of a bulk superconductor vs temperature: T = T c o |l — \ \H\ — Bo\/ H^ 1 } , the dash-dot line 

shows the phase-transition line for a domain-wall superconductivity derived, e.g., in Refs. [H, [(| Here we take B = 1.5H^> . 
(b) Typical plots of superconducting OP distributions presented for different points on the H — T plane: A, B, C and D, where 
the dashed line indicates the domain wall. The vortex chains are shown by the solid circles while the antivortices are presented 
by open circles. 



Substituting the expressions (jT2J) and (fT3")) into the func- 
tional (|4]) and carrying out the minimization over the 
variational parameters 011, 71, a2, 72, and 77 at a fixed 
amplitude Bq of the magnetic stray field we find the crit- 
ical temperature T* 3 vs the applied field H. Typical plot 
of the phase-transition line T* 3 (H) is presented in Fig. [9] 
by a solid line for Bq = 1.5H^ . One can see that apply- 
ing an external magnetic field H we obtain an increase in 
the critical temperature T* 3 of localized superconductiv- 
ity due to the partial magnetic field compensation inside 
the domains. 

As a next step we restrict ourselves to the case H = 
and analyze the dependence of the critical amplitude Bq 
of the domain stray field corresponding to the supercon- 
ductivity nucleation vs the tpo angle. Substituting the 
expressions (|12[) and (fT3")) into the Eq. (0| wc carry out 



the minimization procedure over the variational parame- 
ters ai, 71, «2, 72, and r\. Typical plot of the dependence 
Bq((Pq) shown in Fig.[3]by a dash-dot line (3) is in a good 
agreement with the solid line (2) derived within our nu- 
merical simulations discussed below. 

It is important to note that with the decrease in the 
angle ipo we clearly observe an increase in the field Bq. 
This increase obviously occurs due to a partial shrinking 
of the Cooper pair wave function in analogy to the case 
of a superconducting wedge considered above. For small 
angles ipo — > we find: Bq/H C 2 — > 2.2. In the case 
(fQ = ir/2 for an appropriate gauge choice the problem 
of a half - infinite domain wall can be exactly mapped 
to the wedge problem for x = ft/2. Indeed, for a domain 
wall perpendicular to the edge T the GL functional is 
symmetric with respect to the parity transformation x — ¥ 
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—x and, thus, the OP wave functions are either odd or 
even in the x variable. The even solutions corresponding 
to the energy minimum clearly satisfy the condition 



of the time dependent GL model 



dV(r)/da 







(15) 



which coincides with the one imposed for the supercon- 
ducting wedge at the side with x — 7r /2- As a con- 
sequence, we find B*(n/2) = H™(n/2) ~ 2H c2 . The 
variational parameters corresponding to this particular 
case should be taken as follows: a.\{n/2) = a2(7r/2) = 
0.14, 7x(7r/2) = 72(tt/2) = -0.5, ^(tt/2) = 0.6. 

However this mapping between the domain wall and 
the wedge, of course, does not hold exactly for ipo ^ ir/2. 
Still similarly to the wedge the shrinking of the wave 
function at small angles ipo can result in the appearance 
of a vortex - free solution (see Sec. II C). 

Two domain walls. The above variational results for 
the crossing domain wall and the sample edge remain 
valid also for two domain walls crossing at rather large 
angle ipo < tt/2. Due to the symmetry of the magnetic 
held profile, this generalization is straightforward. 



III. NUMERICAL SIMULATIONS 

To confirm the findings obtained within the trial func- 
tion approach we proceed with the numerical analysis of 
the superconductivity nucleation in a thin-film F/S bi- 
layer within the time - dependent GL formalism. Let 
us consider a superconducting thin-film rectangle in the 
plane (xy) of the cartesian coordinate system (having the 
lateral dimensions L x W) in the presence of a homoge- 
neous external magnetic field H normal to the film plane 
and the step-like magnetic field b(x) = zo-Bo sign(x) of a 
straight domain wall oriented at a certain angle ipo with 
respect to the film edges y = ±W/2 (see Fig. I10p. 




FIG. 10: (Color online) Schematic view of a F/S bilayer: a 
superconducting rectangle (L x W) and a straight domain 
wall (dashed line) oriented at a certain angle <^o with respect 
to the edges y = ±W/2. Bo is a stray field amplitude inside 
the domain, H is an external magnetic field, the x axis is 
chosen to be perpendicular to the domain wall, T is the film 
boundary. 

In order to obtain stationary OP distributions we sim- 
ulate the relaxation to an equilibrium state on the basis 



d 

h— + 2ie${r) 



\&(r) = a*(r) 



f|*(r)|^ (r) + -L 



-ihV - —AM 
c 



divj„(r) 



*(r) , 
divj s (r) , (16) 



where the parameter r\ controls the rate of the OP relax- 
ation, $(r) is the electrochemical potential, a and j3 are 
the parameters of the GL expansion, j„(r) = — <jVQ(r) 
and 



J.(r) = ^ i *(r) 



2e 

-ifiV A(r) 

c 



*(r) 



(17) 



are the densities of normal and superconducting currents, 
a is a normal conductivity, c.c. stands for complex conju- 
gation. Focusing on the study of superconductivity nu- 
cleation one can neglect the effect of vanishing supercur- 
rents on the magnetic field distribution and consider a 
fixed vector potential profile corresponding to the mag- 
netic field B z (r) = H + b z (x). 

The absence of normal current through the sample 
boundary imposes a boundary condition on the poten- 
tial $(r): 



<9$(r) 



On 







where n is a unit vector normal to the edge T of the 
sample. The calculations have been made for two types 
of boundary conditions on the superconducting OP wave 
function: (i) superconductor - insulator boundary con- 
dition 



\-ih^- - -An]*(l 
L on c 1 



, 



(18) 



and (ii) superconductor - normal metal boundary condi- 
tion 



*(r) 







(19) 



To study the equilibrium phase diagram we need all 
transient processes to be finished which corresponds to 
zero electrochemical potential ^(r) = and zero time 
derivatives 94'(r,t)/9t = 0. In our calculations we 
stopped the simulation procedure when the maximum 
of the electrochemical potential reaches the accuracy of 
the numerical computations (10 -15 x ft|e|ao/2mcr/3 in the 
dimension units), where a® — \a(T = 0)|. The calcula- 
tions have been carried out for a grid size 0.2£o x 0.2£o and 
77 = 10X ma/3/he 2 . The time interval between two subse- 
quent iterations was chosen as follows: 0.01 x ma[3/e 2 ao. 

Let us start with the simplest case of a superconduct- 
ing square (L = W = 20£o) placed only in a homoge- 
neous magnetic field H. We consider here the boundary 
conditions for ^ function in the form (|18[) . By varying 



10 




FIG. 11: (Color online) The OP distribution in a superconducting square (20£o x 20£o) in an increasing external magnetic field 
H. The boundary condition (|18[) was imposed on the ^ function. Solid circles show the vortex chains. 




FIG. 12: (Color online) The OP distribution in a superconducting square (20£o x 20£o) for an increasing amplitude Bo of the 
magnetic stray field in the domains and tp — n/2. The boundary condition (|19|) was used for function in order to suppress 
the superconductivity nucleation near the corners of the sample. The vortex and antivortex chains are shown by solid and open 
circles, respectively. 



the external field H we study the evolution of the OP 
distribution ty(x,y) in the superconducting square (see 
Fig. [TTj) . One can see that the Cooper pair wave func- 
tion remains nonzero in the vicinity of the sample corners 
at the field above H c %. Obviously, the corresponding nu- 
clei at the corners can be considered separately only for 
a rather large L values well exceeding the nucleus size. 
In this case each nucleus describes the OP distribution in 
a superconducting wedge with a corner angle ir/2. The 
corresponding nucleation field H™ 3 ~ 2H C 2 appears to be 
in a good agreement with the results of our variational 
analysis in Sec. II D. Our numerical simulations also give 
evidence for the appearance of the vortex chains intro- 
duced in the above consideration. These chains outgo 
from all four wedge vertices along the corresponding bi- 
sectrices (see Fig. [TTj) . 

We continue with the numerical analysis of the super- 
conductivity nucleation in a superconducting rectangle 
L x W affected in the field of a straight domain wall. In 
order to suppress the superconductivity nucleation near 
the sample corners we use mixed boundary conditions for 



the ^ function: the conditions (|18p and (|19p are taken at 
the edges y = ±W/2 and x = ±L/2, respectively. Shown 
in Fig. [12] is the transformation of the OP distribution 
for a square with L = W = 20£o <fo — n/2 caused by 
the increase in the amplitude Bq inside the domain. The 
numerical results are in a full agreement with our ana- 
lytical findings in Sec. II D: (i) the superconductivity is 
localized near the crossing of the domain wall and the 
sample edges y = ±W/2 and survives up to the critical 
value Bq = H™ ~ 2H c2 , (ii) we observe the vortex (an- 
tivortex) chains outgoing from the crossing points of the 
superconducting channels. 

Applying an external magnetic field H we find out 
the phase diagram for localized superconducting states. 
For a particular case Bq = 1.5H^ the corresponding 
phase transition line is shown in the Fig.|H]by the dashed 
line. We also observe the transition lines corresponding 
to the bulk superconductivity (dotted line in Fig. [9]) and 
to the domain wall superconductivity (dot-and-dash line 
in Fig. Typical contour plots of the OP distribu- 
tions for different parts of the phase diagram (see Fig. ^ 
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illustrate the switching between the bulk and localized 
superconductivity nuclcation. 



In the limit of zero external field H we have also ana- 
lyzed the dependence of the critical field amplitude Bq on 
the ipo angle for the sample sizes L = 30£o and W = 15£o- 
Typical plot of the dependence Bq(<Po) is shown in Fig. [3] 
by a solid line (2). Both for small and large corner an- 
gles the numerical dependence is in a good agreement 
with the results of the variational analysis carried out in 
Sec. II D. Typical contour plots of the Cooper pair wave 
function presented in Fig. [T3"l show the transformation of 
the vortex patterns with the changing cpo angle. 



Let's focus on the vortex (antivortex) arrangements 
presented in Figs. [H QTJ [121 H31 Near the crossing points 
of the superconducting channels the vortex patterns ob- 
tained from our numerical simulations appear to be in 
a good agreement with the variational predictions for 
the infinite superconducting channels (see Sec. II) where 
the vortices (antivortices) form the vortex (antivortex) 
chains. Far from these crossing points the vortex distri- 
butions are more complicated and strongly depend on the 
sample geometry and corresponding boundary conditions 
for the Cooper pair wave function. 



IV. CONCLUSION 

To sum up, we have investigated the distinctive fea- 
tures of superconducting OP nucleation for the inter- 
acting superconducting channels in strong magnetic field 
H > H c2 . We have studied three generic problems: (i) 
the OP nucleation between two superconductor/vacuum 
boundaries forming a superconducting wedge; (ii) the OP 
nuclcation between domain wall and the sample edge; 
(iii) the OP nucleation between two domain walls. We 
have shown that in all these cases the crossing of local- 
ized modes results in the increase in the superconducting 
critical temperature. Using both numerical and varia- 
tional analysis of these problems we have developed a 
description of the interference phenomena which govern 
the structure of the OP patterns. The resulting critical 
temperature enhancement and its magnetic field depen- 
dence should be observable in resistive measurements of 
hybrid F/S structures. 
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